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75. 
ON THE ATTRACTION OF AN ELLIPSOID. 


[From the Cambridge and Dublin Mathematical Journal, vol. tv. (1849), pp. 50—65.] 


Part I.—On LEGENDRE’S SOLUTION OF THE PROBLEM OF THE ATTRACTION OF AN 
ELLIPSOID ON AN EXTERNAL POINT. 


I PROPOSE in the following paper to give an outline of Legendre’s investigation 
of the attraction of an ellipsoid upon an exterior point, [“ Mémoire sur les Intégrales 
Doubles,” Paris, Mem. Acad. Sc. for 1788, published 1791, pp. 454—486], one of the 
earliest and (notwithstanding its complexity) most elegant solutions of the problem. 
It will be convenient to begin by considering some of the geometrical properties of a 
system of cones made use of in the investigation. 


§ 1. The equation of the ellipsoid referred to axes parallel to the principal axes, 
and passing through the attracted point, may be written under the form 
l(w@—ayl+m(y—bP+n(z—cP-—k=0, 


} 


(where y, A af A T 4 are the semiaxes, and a, b, c are the coordinates of the 


attracted point referred to the principal axes). Or putting la?+mb?+n?—k=6, this 
equation becomes 
| la? + my? + ne — 2 (lax + mby + ncz) + ô= 0. 

The cones in question are those which have the same axes and directions of 
circular section as the cone having its vertex in the attracted point and circumscribed 
about the ellipsoid. The equation of the system of cones (containing the arbitrary para- 
meter w) is 

(la? + my? + n2?) ò — (law + mby + nez? + œ (2 +y + 2) =0; 
or as it may also be written, 


(@? + 15 — Pa?) x + (w? + md — mb?) y? + (w + nb — n?) 2 — Amnbeyz — 2nleazw — 2lmabay = 0. 
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For w=0, the cone coincides with the circumscribed cone; as œ increases, the 
aperture of the cone gradually diminishes, until for a certain value, œ =Q, the cone 
reduces itself to a straight line (the normal of the confocal ellipsoid through the 
attracted point). It is easily seen that Q? is the positive root of the equation 


Pq? mb? nc? 
0418 pmi Pina” 
a different form of which may be obtained by writing on", E being then determined 
by means of the equation 
oiensireh miè daioun 
k+l ktm kng ~’ 


that is, VG+E) J (+8): J (+8) 


are the semiaxes of the confocal ellipsoid through the attracted point. 


In the case where œw remains indeterminate, it is obvious that the cone intersects 
the ellipsoid in the curve in which the ellipsoid is intersected by a certain hyperboloid 
of revolution of two sheets, having the attracted point for a focus, and the plane of. 
contact of the ellipsoid with the circumscribed cone (that is the polar plane of the 
attracted point) for the corresponding directrix plane: also the excentricity of the hyperboloid 


is Leat mb? + n?e), which suffices for its complete determination. For w=0, the 


hyperboloid reduces itself to the plane of contact of the ellipsoid with the circum- 
scribed cone, and for œ =Q, the hyperboloid and the ellipsoid have a double contact, 
viz. at the points where the ellipsoid is intersected by the normal to-the confocal 
ellipsoid through the attracted point. 


If œ remains constant while Æ is supposed to vary, that is, if the ellipsoid vary in 
magnitude (the position and proportion of its axes remaining unaltered), the locus of 
the intersection of the cone and the ellipsoid is a surface of the fourth order defined 
by the equation 

(la? + my? + ne — lax — mby — nez} = œ (a? + yY? +2), 


and consisting of an exterior and an interior sheet meeting at the attracted point, 
which is a conical point on the surface, viz. a point where the tangent plane is 


replaced by a tangent cone. The general form of this surface is easily seen from the 
figure, in which the ellipsoid has been replaced by a sphere, and the surface in question 
-© 55 
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is that generated by the revolution of the curve round the line CM. The surface of 
the fourth order being once described for any particular value of œ, the cone corre- 
sponding to any one of the series of similar, similarly situated, and concentric ellipsoids 
is at once determined by means of the intersection of the ellipsoid in question with 
the surface of the fourth order. It is clear too that there is always one of these 
ellipsoids which has a double contact with the surface of the fourth order, viz. at the 
points where this ellipsoid is intersected by the normal to the confocal ellipsoid 
through the attracted point; thus there is always an ellipsoid for which the cone 
corresponding to a given value of w reduces itself to a straight line. 


Consider the attracting ellipsoid, which for distinction may be termed the ellipsoid 
S, and the two cones C, C’, which correspond to the values œ, »—dw of the variable 
parameter. Legendre shows that the attraction of the portion of the ellipsoid S 
included between the two cones C, C” is independent of the quantity k, which determines 
the magnitude of the ellipsoid: that is, if there be any other ellipsoid 7’ similarly situated 
and concentric to and with the ellipsoid S, and two cones D, D’, which for the 
ellipsoid T correspond to the same values w, w— dw of the variable parameter; then the 
attraction of the portion of the ellipsoid S, included between the two cones C, C’, is 
equal to the attraction of the portion of the ellipsoid T included between the two 
cones D and D’. By taking for the ellipsoid 7 the ellipsoid for which the cone D 
reduces itself to a straight line, the aperture of the cone D’ is indefinitely small, and 
the attraction of the portion of the ellipsoid T included within the cone D’ is at once 
determined ; and thus the attraction of the portion of the ellipsoid S included between 
the cones 0, C’ is obtained in a finite form. Hence the attraction of the portion 
of the ellipsoid S included between any two cones C,, C, corresponding to the values 
w, œw, of the variable parameter, is expressed by means of a single integral, and by 


extending the integration from w=0 to w= T E) the attraction of the whole 


ellipsoid is obtained in the form of a single integral readily reducible to that given 
by the ordinary solutions. It is clear too that the attraction of the portion of the 
ellipsoid S included between any two cones C,, C,,, is equal to that of the portion 
of the ellipsoid T included between the corresponding cones D, and D,. Hence also, 
assuming for the ellipsoid T, that for which the cone D, reduces itself to a straight 
line, and supposing that the cones C, and D, coincide with the circumscribing cones, 
the attraction of the portion of the ellipsoid S exterior to the cone C, is equal to 
the attraction of the entire ellipsoid T. More generally, the attraction of the 
portion of the ellipsoid S included between the cones C, and C, is equal to the 
attraction of the shell included between the surfaces of the two ellipsoids, for which 
the cones D, and D, respectively reduce themselves to straight lines. 


§ 2. Proceeding to the analytical solution, and resuming the equation of the 
ellipsoid 
læ + my? + n2* — 2 (lax + mby + nez) +8 =0, 


and that of the cone 


(la? + my? + ne?) ò — (lax + mby + nez} + œ (2 +y + 2) =0; 
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consider a radius vector on the conical surface such that the cosines of its incli- 
nations to the axes are 

ne a 

© > e > O’ 
P,Q, R and © being functions of the parameter œ, and of another variable œ, which 
determines the position of the radius vector upon the conical surface. Also let p be 
the length of the portion of the radius vector which lies within the ellipsoid; then 
representing by dS the spherical angle corresponding to the variations of œ and ¢, the 
attraction in the direction of the axis of æ is given by the formula 


P ay 
Also by a known formula 


-aP oe) 0(F ete as) Beas sat 


O=V(P+ @ + Rj, 


dp do ôw dd 


dp dw dd dw 
and it is easy to obtain 
k? 200? 
P = 1P? + mEt nk’ 


The quantities P, Q, R have now to be expressed as functions of œ, ¢, so that 
their values substituted for æ, y, z, may satisfy identically the equation of the cone. 
This may be done by assuming 


P=, 
Q=mb (w° + nò) (+5 +, cos )-2 ZNP ne sin d, 
R=ne (œ + mò) (la+ +5, f cos $) += oe mb sin ¢, 


where 
P = (w + Mò) (w + nd) — mh? (w + nb) — n? (w? + mS), 


U? = mh? (w? + nd) + ne (w? + mò), 


mb? nè 1 


D? = (a? + 18) (w? + mò) (o +23) oes aple att md* 4nd j 


a system of values which, in point of fact, depend upon the following geometrical con- 
siderations: by treating æ as a constant in the equation of the cone, that is, in effect 
by considering the sections of the cone by planes parallel to that of yz, the equation 
of the cone becomes that of an ellipse; transforming first to a set of axes through 
the centre and then to a set of conjugate axes, one of which passes through the 
point where the plane of the ellipse is intersected by the axis of æ, then the equation 
2 2 
takes the form ia and is satisfied by E= A cos ġ, n= Bsin ġ, and z s : being 
of course linear functions of these values, the preceding expressions may be obtained. 
55—72 
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The substitution of the above values of P, Q, R (a somewhat tedious one which 
does not occur in the process actually made use of by Legendre) gives the very 
simple result, 

_ Pio dodo, 
dS = Wig: 


and the formula for the attraction becomes 


ye mn WO a Pto do db 
IP + m@ + nR? 


which is of the form A =2 | Iw*dw, where 


Pei Phd¢ 
=] Pime + 0R” 


which last integral, taken between the limits ¢=0 and ¢=27, and multiplied by 
2w*dw, expresses the attraction of the portion of the ellipsoid included between two 
consecutive cones. The integration is evidently possible, but the actual performance of 
it is the great difficulty of Legendre’s process. The result, as before mentioned, is 
independent of the quantity k, or, what comes to the same thing, of the quantity ò: 
assuming this property (an assumption which in fact resolves itself into the consideration 
of the ellipsoid for which the cone reduces itself to a straight line, as before explained), 
the integral is at once obtained by writing 6=A where A represents the positive root 
of the equation 
Pa? mh? nc? l 0 
w+ ote mA  atend oT 
This gives 
_ Pa? (œ + mA) (wo? + nA) 
w + lA $ 
Q = lmab (œ + nA), 
R=Inac (œ+ mA), 


} 


values independent of ¢, or the value of J is found by multiplying the quantity under 
the integral sign by 27: and hence we have 


w (a? + IA)? (o + mA)! (a? + nA)! do 


A = dela Car mB E MAF Cot 4 nA CoP FID 00 (OFF ROY (OF FTO)” 


where of course A is to be considered as a function of œ. By integrating from 
w=, to o=@,, we have the attraction of the portion of the ellipsoid included 
between any two of the series of cones, and to obtain the attraction of the whole 
ellipsoid we must integrate from w=0 to œ = Hi @) , where & is determined as before 
by the equation 

la? mb? nc 


bie sae? CA 
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and it is obvious that for this value of œ we have A=6 The expression for 


the attraction is easily reduced to a known form by writing y ="; this gives 


A= 4rla | K (k + ly} (k + my)! (k + nyodo ; 
Ba? (le + my) (le + ny) + mb (le + ny) (e +y) + niet Q + ly) (e+ my 


Also i 
ee k( Pa? = mb? n nec b 
— \k+ly k+my k+ny/’ 
whence 
EEN k [Ba? (k + my)? (k + ny)? + mb (k + ny)? (k +y} + nic? (k + ly) (k + myy] 
2 (k + lyf? (k + myY (k+ nyy 
and thus | A=2ark'la I dy ’ 
mt (k + ly)! (k + my)? (k + ny)? 


where for the entire ellipsoid the integral is to be taken from y=& to y=o. A 
k+lé 


, in which case the limits 
k+ly 


better known form is readily obtained by writing a= 


for the entire ellipsoid are æ= 0, x=1. 


It may be as well to indicate the first step of the reduction of the integral J, 
viz. the method of resolving the denominator into two factors. We have identically, 


(A — 8) (LP? + mÆ +n R?) = œ (P?+ Q?+ R?) + A (LP? + mÆ + nE) — (laP + mbQ + ncRy, 


and the second side of this equation is resolvable into two factors independently of 
the particular values of P, Q, R. Representing this second side for a moment in the 
notation of a general quadratic function, or under the form 


AP? + BÆ + CR? + 2FQR + 2GRP + 2HPQ, 
we have the required solution, 


LP? + mQ? + nR? = 
JAP + {H+ V(- ©) Q + (4+ v- 38)} R][AP + {H - v(— ©) Q + (G — v(— 38)} R]; 


where, as usual, 3 = CA - @, @ = AB- H?, and the roots must be so taken that 
V-B) J- ©) = f (F = (GH - AP). 
I have purposely restricted myself so far to the problem considered by Legendre : 


the general transformation, of which the preceding is a particular case, and also a 
simpler mode of effecting the integration, are given in the next part of this paper. 
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Part I]—On A FORMULA FOR THE TRANSFORMATION OF CERTAIN MULTIPLE INTEGRALS. 


Consider the integral 
v=| Fo y, ...) de dy ..., 
where the number of variables æ, y,... is equal to n, and F(s, y,...) is a homogeneous 
function of the order p. 


Suppose that a, y,... are connected by a homogeneous equation (a, y, ...)=0 
containing a variable parameter œ (so that w is a homogeneous function of the order 
zero in the variables a, y,...). Then, writing 


P=H+yt+..., @=7a, y=rP,... 
the quantities a, 8, ... are connected by the equations 
COLO ES.= Li a O 8... y= 6, 


and we may therefore consider them as functions of œw and of (n—2) independent 


variables 0, ¢, &c.; whence 
da dy...=1""V drdwdé..., 


where 
Taia Bop a 
da dp 
dw’ da’ 
da dg 
dé ? dé’ 
Also F(a, Y.. JEG P- 


and therefore 
V= |m F(a, B,...)V dr dwd@ ..., 
or, integrating with respect to 7, 


| petn- dr a 1 gate. 
pen 


which, taken between the proper limits, is a function of a, £,..., equal f(a, £, Pos) 
suppose; this gives 


V=|fla B, ...) F(a, B, ...) V do dð... 


in which I shall assume that the limits of œ are constant. If, in order to get rid 
of the condition a+ 8? ...=1, we assume 


p q 
a=}, B=4,..., p=pt¢t... 
p $ p= pi +q 
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the preceding expression for V becomes 


v=|r(. 1, ...) Fp, DE Sa 


in which 
AEE ET A WENN vials 
dp - dq 
dw’ dw’ 
dp dg 
dé’ dé’. 
Assume 


p=PE+ Pnt+P%... 
g =QE +n +E.. 


where the number of variables £, 7, €... (functions in general of œ, 0, &c.) is n, and 
where the coefficients P, Q, &c. are supposed to be functions of œ only. We have 


dp _ dé , dn , dE 
da= t got? Pasa 


e ET ARN Oo 


and, substituting these values as well as those of p, q, &c., but retaining the terms 


e L, &c. in their original form, the determinant D resolves itself into the sum 


of a series of products, 


| d dg lilura a6 on 
dw’ dw’ 

A ii eae! ae 
Pe , Q ? dn dt 
Pr. Oi. . de’ dé’ 


Let W be the function to which w (p, q, ...) is changed by the substitution of the 
above values of p, q,... so that Y is a homogeneous function of £E, n, ¢,... and we have 
the relation Y=0. (It will be convenient to consider é, n, &,... as functions of 
w, 0, &e., such as to satisfy identically this last equation.) We deduce 


dé’ dé’ dé’ ° dé’ 


res Pe E AY lt 1, . ... |= &a = 8 suppose, 
a d n ? gc ? Y E ? bi £ ? 
| dn dý dé dg 
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where for shortness X = oe a A , &c. The substitution- of these values gives 
D=| dp dg S 
dw’ dw’ |’ 
Aion ye | 
fo. a Se 
ay Mew Beas 


where it will be remarked that the successive horizontal lines (after the first) of the 
determinant are the differential coefficients of Y, p, q,... with respect to & with respect 
to n, &c. In general, if y denote any function of p, qg,... these quantities being 
themselves functions of œ, £E, 7,..., and £, n,... containing w; also if W be what y 
_ becomes when for p, g,... we substitute their values in a, &, n, ...; then we have identically 


is Pa 


dw’ do’ oe =0. () 
AS E A 
Toa W 

Q”, 


Z WR 


In the present case however, writing for shortness W(p,q,...)=~, this function y 
contains w explicitly as well as implicitly through p, q, &c. The formula is still true if for 


oy we substitute nit il o a ; 7 on the second side denoting a partial differential 
dw dw dw’ dw 


coefficient taken only so far as œ is explicitly contained in y. And considering p, g,... as 
functions of œ, E, 7,... (E, 7,... themselves functions of w and of other variables which 


need not here be considered), y or VY vanishes identically, and we have a 0. Hence, 


dy dy 


in the last formula, we have to write wai instead of Hos’ and we thus derive 
dy|P, Q...|=| dp dq |, 
dw P Q, . dw ? dw pee 
F X , ie , Q ) 
Vy Fo aS 


Ai, Po, gs 


1 This formula, or one equivalent to it, is given in Jacobi’s memoir ‘‘De Determinantibus Functionalibus,” 
Crelle, t. xxi. [1841] p. 319. 
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whence D becomes 7 
a ave 
D=| P, Q, | FES, 
P, Q, 


which is the value to be made use of in the equation 


v-|7(, 1...) Fe, dr) sen D ded... 


The principal use of the formula is where y is a homogeneous function of the second 
order of p, g,.... Thus, suppose 


y= 4 (4p + Be... +2Hpqt+...), 


also, for the sake of conformity to the usual notation in the theory of transformation 
of quadratic functions, writing a, a’,... B, B’, ... instead of P, P’,...Q, Q,... and putting 
after the differentiations £=1, we have 


p= a+ an +a” E+ oeo j 
q=B+Bn+ REH., 


values which we may assume to give rise to the equation 
(dœ + Be... +2Hpq ...)= (1 -P &-...), 


(where 7, & ... are taken to be functions of 0, &c. such as to satisfy identically the 
equation 1 = 77+ f?+...). . 


Hence, by a well-known property, if 


we have 
Ce Pee = JE 
th r 

so that, observing that in the present case X = 1, and therefore 


BE a E 


dn dt 
dé’ d@ 
we have 
heer i ale 
rafe re 2, ) FO, D o) Sern To S U 
C 56 
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The remainder of the process of integration may in many cases be effected by the 
method made use of by Jacobi in the memoir “De binis quibuslibet functionibus &c.” 
Crelle, t. xu. [1834] p. 1, viz. the coefficients a, a’, &., B, &e., may in addition to the 
conditions which they are already supposed to satisfy, be so determined as to reduce 
any homogeneous function of p, q, 7,... entering into the integral to a form containing 
the squares only of the variables. This method is applied in the memoir in question 
to the integrals of n variables, analogous to those which give the attraction of an 
ellipsoid; and that directly without effecting an integration with respect to the radius 
vector. I proceed to show how the preceding investigations lead to Legendre’s integral, 
and how the method in question effects with the utmost simplicity the integration 
which Legendre accomplished by means of what Poisson has spoken of as inextricable 
calculations. 


Consider in particular the formula 


yaf 2 ii) ADAY: sss 
Cty. P oi 


the number of variables being as before n, and (s, y,...)* denoting a homogeneous 
function of the order h. The equation for the limits is assumed to be 


l(æ—a}+m(y—b?} +... =k. 


Assume 
ba, Y, ...)= 0? (+ y+...) + Ò (le? + my + ...)— (lax + mby + ...}, 


(where ô, = la? + mb... — k, is taken to be positive); or more simply, 
y (æ, y...) = (w + 18 — la’) a + (w+ mò — mb) y? +... — Alm ab sy —... 
Here »=h+2i—3, Also, putting for shortness 
lap + mbq +... =A, lp?+mq+...=®, 


it is easy to obtain 


cp feeds, Ae bovo $ nting p ANE E etl ee Ser ra 
p’ p ~~ ht+24+n—8 pitzitn—s 
TSA: d 


Pa? mb? 


r = (w? + 18) (wo? + md)... C 


5 &e.) 


values which give 


{ise —)tin—) wp" [( A+ wp) titn—s aa ( AS p)rtaitn—s] (p, Oia i Sdo dé ... 


Dprtzi+n—s 


= 3° 
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where it will be remembered that p, q,... are linear functions (with constant terms) 
of (n—1) variables n, ¢,..., these last mentioned quantities being themselves functions 
of (n—2) variables 6, &c. such that 1—7?-&-.. ai aon If besides we suppose 


that = lp? + mg? + ... reduces itself to the form r — 5- &c., we have, by the formula 
of the paper “On the Simultaneous Transformation of two Homogeneous Equations of 
the Second Order,” [74], 


(30-3 
1 
Plamen 18 — Dr) (w + mò — md) ons (1 “a ba emp mse mn = se 
which is true, whatever be the value of X. 
It seems difficult to proceed further with the general formula, and I shall suppose 
n=8, t=0, h=1, (x, y...)*=a, or write 
y= [See æ da dy dz 
CEPT 
the equation of the limits being 
l (x —a}+m(y—bf +n (z— c} =k. 
Here we may assume ņ=cos0, ¢= sin 0, (values which give S=1). And we have 


V=2 wdw Seah a cos 0 + a” sin 0) do. 
5 cos? 0 sin? 0 


-— 


z+ Q R 


from 6=0 to 0=2r; or, what comes to the same thing, 


w dw ad@ 
Lis Je |1 asoa co 0 sin?’ 


gH, HEE 


ye y 
from 0 =0 to 0=4r. Hence 


awPy(QR)do 
VP -Q P-B) 


we have from the formulæ of the paper before quoted, 


QR (B- mP) (C -nP)- P 
K (P-Q) (P-R) 


B, O, F, being the coefficients: of y*, 2°, yz in ẹ (e, y, z), viz. 
B= +m- mb, C= +n- re, F=mnbe; 


v=4r | - 


a = 


and consequently PQR {(B—mP) (C -— nP) - F°} 
V =4r | o do = ra (P-Q (P-R) 


56—2 
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Also from the equation 


(3) 6-9)0-2)- 


} (w? + 18 — 1d) (œ? + mò — mA) (w? + nS — nr) (1 Pa? mb? nec? ) 


wte w+md—mr FNA 
differentiating with respect to à, and writing A =P, 
1 
-p(t- 9) (1-7) 
(œ + 15 — LP) (w° ++ mò — mP) (w + nd — nP) 
Ae mb? nic? ) 


1 
=~ wt lò Py (co? + md — mPF (o + nd —nP){’ 
or, as this may be written, 


Bae 5. | See 
«(P—Q)(P-B) 
lh 1 
(a+ WIP) (ot mB —mP) (ot nd —nP) Carp + Core mP Ft FREE 


and from the values first written down, for B, C, F, we obtain (B—mP) (C—nP)—F 


=(w?+md)(w?+n8) —mb?(w? +d) —n*e?(w? + md) -mP (w* + nbd —n*c?) —nP (wm — mb) + mnP* 
=(w? + md — mP) (w? + nd — nP) — mb (w + nb — nP) — nic? (w° + md — mP) 


= Dal Co tee ee Ane e iad, the last reduction being effected by means of the equation 


w° + lò — LP i 
1 ee ee T ee a 
w + lè -lP wt+md—mP o?+nd—nP ~*~ 
Hence 
PQR {(B — mP) (0 — nP) — F3} 
PA T A p: 
la 


 (w*+18—1P)! mP pioppi ed Bes off Ao atlas ay ole 
(w?+16 LP) (w?-+-mé mPY (w?+nd—nP) VIP + lemi —=mPP a (o*-+nd—nP}} 
Substituting this value, and multiplying out the fractions in the denominator, 


V = 4rla x 
| w (œ? + 18 — LP) (w? + mS — mPÈ (w? + nd — nP} do | 
Ba? (@?-+md—mP)(o?+nd—nP)? + mb? (w?+-nd—nP)? (o+ -IPY +e (o+ -LPY (w+ md— mP} 


the reduction of which integral has been already treated of in the former part of 
this present memoir. 
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